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ABSTRACT 


In Part A, systematic nunerical solutions of two dimensional 
and axisymmetrical laminar jet of an incompressible fluid with and 
without free stream have been obtained. The exact numerical solutions 
have been checked with experimental results and similarity solutions 
for the case without free stream. At far downstream, the numerical 
solutions approach the values of similarity solution. The numerical 
solutions give better agreement with experimental data than the simila- 
rity solutions. In general, it shows that the boundary layer equation 
is a good approximation of laminar jet problem provided that the Reynolds 
number at nozzle exit is not too low. With free stream, the numerical 
solutions agree with linearized analytic solution if the jet excess 
velocity is small in comparison with the free stream velocity. The. 
non-linear effects are to decrease the rate of decrease of central 
velocity of the jet and to broaden the spread of the jet. 

In Part B, an approximate numerical solution for three dimen- 
sional laminar jet is proposed. The accuracy of the method has been 
determined from the exact solutions of the two limiting cases of three 
dimensional jets, i.e., two dimensional and axisymmetric case of Part A. 
It was found that this approximate method may give good results for the 
axial velocity distribution and the spread of three dimensional jet. 


SYMBOLS 


Radius of circular exit 


Half-width of jet with and without cross flow respectively 

Momentum coefficient for two-dimensional and axlsymmetric jet 
respectively 

Correction factor for central excess velocity 
Half-width of the shorter side of rectangular jet exit 
Spacing in x direction 

Bessel function of zeroth and first order of the first kind 
respectively 

Spacing in z direction 

Aspect ratio of rectangular jet exit 

Total momentum of two dimensional and axlsymmetric jet respec- 
tively 


Spacing in y direction * 

U.D 


Reynolds number defined 


j 


V 


and 


U a 


respectively 


* * * 

Excess velocity components along x , y and z direction respec 


tively 


* . * 


Non-dimensional excess velocity components defined by u /U^ , 


R v /u . and R w /U. respectively 
e j e j 



%v 


SYMBOLS 




U 1 

A 

u 


u 

u 


o 

A 

j 



AAA 

x , y , z 


x, y, z 


5 

6 

v 

P 


Non-dimensional maximum excess velocity over a station with 
and without cross flow respectively 

Non-dimensional excess velocity distribution at jet exit 


Uniform velocity of moving stream 

•fa ^ 

Non-dimensional uniform velocity of moving stream (U /U.) 

o' j 


Excess velocity at the center of jet exit 


Mean excess velocity at jet exit 
Axis of Cartesian coordinate 

Non-dimensional distance defined by x /R D , y /D and 
a e 

z /D respectively 

Non-dimensional parameter defined by Eq. (10) 

A constant 

Kinematic viscosity of fluid 
Density of fluid 


Subscript: 

o Axisymmetrical jet 

o° Two-dimensional jet 

L Three dimensional jet with aspect ratio L 

Superscript : 


Without cross flow 


I. INTRODUCTION 


In most of the theoretical Investigations of laminar jet 

1 2 

mixing, only the similar solutions 9 or the solutions of linearized 

1 3 

equations * have been discussed. Very little general solutions of 
laminar jet mixing have been obtained. In the first part of this 
report, we study numerical solutions of two dimensional as wall as 
axisymmetric laminar jet mixing of an incompressible fluid with and 
without free stream sys temmaticaliy , We are especially interested In 
the non-similar solutions of these problems. The Influence of initial 
velocity profiles and the effect of non-linearity will be analyzed, , 

The results of our numerical solutions will be compared with the well 
known similar solution (without free stream) and linearized solution (with 
free stream"*). 

In the second part of this report, numerical solutions of a 
three dimensional laminar jet mixing have been studied. In our preliminary 
attempt of solving the complete non-linear equation of a three dimensional 
jet mixing by computer, the numerical calculation was very unstable and so 
far no satisfactory results were obtained. Since we have satisfactory so- 
lutions for two dimensional and axisymmetrical laminar jet mixing in Part I, 
we try to find some approximation in numerical solution in the three dimen- 
sional case from the results of our two dimensional and axisymmetrical cases. 
We find that the two dimensional (aspect ratio is infinite) and the axi- 
symmetrical case (aspect ratio is unity) represent two limiting cases of 
three dimensional laminar jet mixing from a rectangular nozzle. Hence in 


// 


2 


general, the results of three dimensional jet mixing lie between those 
of two dimensional and axlsymmetric cases. From the results of two 
dimensional and axlsymmetric jets, we found that the axial velocity 
distributions may be obtained with sufficient accuracy by neglecting 
the cross -velocity effects. Hence we use similar approximation by 
neglecting the cross-velocity In solving the axial velocity distribution 
of three dimensional jet mixing, which is the most important quantity In 
practical application of jet mixing problems. We also discuss the accu- 
racy of this approximate method and the improvement pf the results from 
our experience of the exact results of two dimensional and axlsymmetric 
cases, Numerical solutions of axial component of three dimensional jet 
without free stream from rectangular nozzle of aspect ratios 1, 2 and 4 
have been obtained. 
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PART A. TWO DIMENSIONAL AND AX1 SYMMETRICAL JET 


II. FUNDAMENTAL EQUATIONS 

The fundamental equations for a steady two dimensional or an 

t 

axlsymmetrlcal laminar jet mixing of an Incompressible fluid are^’ 
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it 
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< 1 ) 



i *s 

T U ) 




V a 
*6 3y 

y 


(y 


*<\ 3u 


9y 


where u and v are respectively the excess axial 

over the free stream U and the transverse or radial 

o 

component; 6 ** 0 for two dimensional case and 6 » 1 


•) ( 2 ) 

(x* - wise) velocity 
(y* - wise) velocity 
for axlsymmetrlcal 


* 


case; U q Is a constant velocity of free stream and v Is the coefficient 
of kinematlcal viscosity of the fluid which Is assumed to be constant in this 
report. The pressure In the flow field Is assumed to be constant In this 

it 

report. When the surrounding stream is at rest, U Q is zero. 

The Initial and boundary conditions of both two dimensional and 
axlsymmetrlcal laminar jet are; 


* « 
x « 0 ; 

* 

x > 0 ; 


* 

u * 


* 

v *• 


* 

u * 


WjW.y ) 
0 at 

0 at 


* 

y 
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* 

y oo 


( 3 ) 
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We solve Eqa, (1) and (2) with the Initial and boundary 

A ft 

conditions (3) where u^(0,y ) is a given function of y . in the 
numerical calculation, it is convenient to Introduce the following non 
dimensional quantities. 


ft 


u « 
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v r V 7 
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o 
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1 
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D R 
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ft 
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ft 

U .D 


R * 
e v 


(4) 


^ A 

where U, is the value of u at the center line of the exit of the nozzle, 
J 

D is the half width or the radius of the jet exit and R is the Reynolds 

w 

number of the present problem. In terms of the non-dimensional quantities 
of Eq. (4), the fundamental equations (1) and (2) are Independent of the 
Reynolds mffflber R^ and are as follows i 


9u , 9v , » v 
3x 3y y 


(5) 


/T , , V 3u , 3v 19 ,6 3Uv 
(B o + u) 8? + 37 ■ “6 W (y W } 

y 


( 6 ) 


1 9u 3 u 

It should be noted that at the axis of the axisymmetrical case - j 

y y 3y 


and 


9v 

9y 


5 


We are going to solve Eqs. (5) and (6) for the following three 
different initial velocity profiles; 

(a) Rectangular Profile: u^ » 1 y * 1 

(7) 

* 0 y > 1 

(b) Parabolic Profile; u,j ■ (1-y^) y ^ 1 

( 3 ) 

« 0 y > 0 

(c) Triangular Profile: * 1-y y £ 1 

(9) 

u^ * 0 y > 0 


III. SIMILARITY SOLUTIONS OF TWO-DIMENSIONAL AND AXIS YMMETRI CAL JET 
Similarity solutions occur for the case without free stream 

* 

(U ® 0) only. For two dimensional jet, Blckley gave the closed form simi- 

1 

larlty solutions as follows: 

M^ 

u* - 0.4543(-5~) 1/3 tanh % 

p VX 


MnV i y o o 

0.5503(-^) i/J [2$(l-tanIT5) - tanh 

px 


( 10 ) 


^21/3 y* 
g - 0.2752(— i) 1/J -J 


PV 2 ' x * 2/3 


where M 0 is the total momentum across a section of the two dimensional 
jet. 
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For axlsymmetr;lcal jet, Schllchtlng gave the closed foiro 
1 

similarity solution as: 


M 

* Jjr 1_ _ 

“ “ 8 " pvx* (I + 1/4?) 2 



(id 


where H r is the momentum across a section of the axisynimetrical jet. 


Let 

M 2 « C 2 pUj D (12) 

16 

then 53 2, an ^ 2/3 respectively for rectangular, parabolic and 
triangular initial velocity distributions. Similarly, let 

M « C pttD 2 U* 2 (13) 

r r j 

then C “ 1 and 1/3 respectively for rectangular and parabolic initial 
velocity distribution. 

Hen cp, Eqs. (10) and (11) may be written, in terms of non- 
dimensional forms , as follows : 


7 


For two dimensional jet 


C 2 

u - 0.4543(-|) 1/3 (l - tanh (■) 


v - 0.5503(-|) 1 ' ,3 [25(l-tanh 2 5) - tanh 5] 

X 


(14) 


5 - 0.2752 C^ 3 y/x 2/3 . 


For axi-'symmetrical jet 


n C 

3 r 
u B ~ 

8 x 


v ~ 


/3 C 


1/2 


(1 + f« 2 ) 2 


t Ir 3 


4 X ( 1 + i ? 2 ) 2 


A- C 1/2 v 
"¥16 r y 


/x 


(15) 


IV. LINEARIZED SOLUTIONS WITH UNIFORM FREE STREAM VELOCITY U . 

o 

A A 

When U q >> Uj , the fundamental equations of laminar jet 

1 

mixing may be linearized. Pai gave the following linearized solutions; 
For two dimensional jet 


u 


i[ er f (~i^2. ) 

2 2^r 

o 


+ erf( 


1+2L 

2/x 


•)3 


C16) 
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and for axisymmetrlcal jet 


u 


• -X x 

e 0 J (Xy) J (X) dX 

0 


( 17 ) 


where x q * x/U Q . Eqs. (16) and (17) are for the rectangular Initial 
velocity profiles only. 


V. NUMERICAL SOLUTIONS Of TWO DIMENSIONAL AND AXI SYMMETRICAL 
JETS WITH AND WITHOUT UNIFORM FREE STREAM VELOCITY. 

The numerical calculation has been carried out on the Univac 1108 
computer of the Computer Science Center of the University of Maryland. An 
implicit finite difference method as described in Reference 4 was used for 
the present calculation. We write 


ja 

9x 


X - Vi> 


forward difference 


It ■ !? (< Wi ■ w centrai difference 


^2 ” VVl ~ ^m + ^ - i ) central difference 

P 


where n refers to station in x-direction and m * that in y-directionj Q 
denotes any quantity under operation, h = Ax and p * Ay denote respectively 
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the spacing In x and y directions. Substituting the finite difference 
representations in Eqs. (5) and (6), we have 


{1 + [v ' u n,m-l " < 2 + T 5 u n,m + 


+ {1 - [v 


]£} U 


u u 


iv-l»m ^ 


(m-l)p J 2 n,nrt*l 


(13) 


and 


¥ 


6 j 

y r 


O / * 


(19) 


0 


where u and v appear in the coefficient are considered to be the 

known values at station (n,m) . The trapaizoidal rule was used for the 

integration of v . In practical calculation, the previous known values 

of u - and v - were first used in the coefficients for u and 

n-l,m n-l,m 

v to obtain u and v . Then Iteration was carried out by re~ 

n,m n,m ^ 

placing the new u and v in the coefficient until a satisfactory 
r ° n,m n,m J 

result for u and v was obtained. Therefore, the u and v thus 
n,ra n,m 

obtained in a station (n,ra) satisfy the fundamental equations and boundary 
conditions exactly. The method of solving the simultaneous equations for 
U n m ®^ ven Appendix I. 

In principle, the boundary of the free jet is infinite. However, 
in the actual computation, it was found that the boundary could be set at 
20D where D is the half width of the nozzle opening for good results in 
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the range presented In this report. Xn the initial stations, the value 
of h was set at 0,001 and as the computation proceeds downstream, 
the value of h is increased. The values of h used in this report 
are between 0.001 to 0.1 for two-dimensional jet and between 0,001 
to 0.05 for axlsymmetrical jet. Lateral spacing p was always kept 
at a value of 0,005 . A test for p * 0.1 does not change the results 
significantly. 

VI. RESULTS AND DISCUSSION 

In jet mixing problems, we are Interested in the decrease of 
maximum axial velocity downstream, the velocity profiles of both axial u 
and transverse v velocity component and the half jet width which shows 
the spread of the jet. The numerical results for these quantities are 
given below: 

(a) Variation of maximum axial velocity. The non-dimensional 

A 

maximum axial velocities, u » u /U, which occurs on the jet axis, 

m m j 

are plotted as a function of non-dimensional axial distance x with 

rectangular initial velocity profile with various free stream velocity u q 

for two dimensional jet in Fig. 1; while the correspondong curves for 

axlsymmetrical jet in Fig. 2. It should be noted that the non-dimensional 

distances for cases with (x ) and without (x) moving stream are different. 

o 

For the cases with free stream, (U 0 “ ^ ) three cases for U Q - 1/2, 1 and 
10 were calculated. It is Interesting to see that the numerical solution 
approaches the linearized solution (16) or (17) as U q increases. This 
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justifies both the correctness of the linearized theory and the numerical 
solution, 

When U Q « 0 , the numerical solution should be compared with 
the similarity solution (14) or (15) , These comparisons are shown in 
Figures 3 and 4. In Figure 3, two dimensional jets with three different 
initial velocity profiles (rectangular, parabolic and triangular) are 
compared with corresponding similarity solutions (C r « 2, 16/15 and 2/3), 
In Figure 4, axisymmetric jets with two different initial velocity profiles 

(rectangular and parabolic) are compared with corresponding similarity 

» 

solution (C r 83 1 and 1/3), From both Figures 3 and 4, it is seen that 
in each case the similarity solution approaches the corresponding numerical 
solution at far downstream as expected since it is well known that simila- 
rity solution is only good at far downstream. 

There are a few experimental results available for laminar two 
5 6 7 

dimensional * and axisymmetrical jet,' These experimental data are also 
plotted in Figures 3 and 4. 

For two dimensional laminar jet, Reference 5 gives experimental 

data for four different Reynolds numbers. The Initial velocity profiles 

correspond to the parabolic distribution because it is supposed to be 

fully developed velocity profile from a two dimensional channel flow. We 

found that the agreement of experimental data with curve of parabolic 

— * 

initial velocity distribution is fairly well for R £ - 4DU^/3v - 68 and 
31 while for R ■ 240, the data approaches the values with rectangular 
initial velocity profile and for R g * 15, the data points are too low. 
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Since no Initial profiles were given, we suspected that for * 240, 

the Initial profile in the experiment might not reach the fully developed 

stage but still close to entrance rectangular profile. For R *15, the 

velocity of the jet is so low that boundary layer approximations might 

not be valid. From the data R ■ 31 and 68 as well as the data of 

e 

reference 6, R £ « 31.5, the correlation between experimental and nu- 
merical results are excellent. Since in the experimental data momentum 
lost across the section were shown, we have to make the proper correction 
in comparison of the experimental data to numerical solution as follows: 


* 

* 

ZjO 


M, 

d 

M, 


( 20 ) 


20 


where subscript 0 denotes the original value at exit and Mj, is the 

* 

measured momentum at downstream and is the corrected value at exit 

to be used for downstream station. 

In Figure 4, the experimental data of Reference 7 indicate a 

it 

shift of origin of the axis ymmetri cal jet of x g 58 x s /R e D = 0.225 for 

rectangular initial velocity distribution (a maximum value from experimental 

data) and x « x /R D ** 0.1 for parabolic initial velocity distribution 

s s e ft 

(an average value of experimental data. Note in Reference 7, x « 0 . 15 C~~) 
* 

= 0 . 1 (—** — ) u ’ •* -r* u.) . With the shift of origin, it is seen that the 
v j 

experimental data agrees extremely well with numerical solution for parabolic 

initial velocity distribution while it is fairly well for rectangular initial 

distribution. According to numerical result, x for the later should be 

s 


0.3 to give best agreement. 
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From the comparison with experimental data, it shows the 
validity of the boundary layer approximation for laminar two dimensional 
and axlsymmetrica! jet when the Reynolds number is not too low. Numerical 
solution provides the overall information of the whole flow field includ- 
ing those near the exit of the nozzle while similarity solution is good 
only far downstream. 

(b) Axial velocity profile u , The axial velocities u vs. 
the transverse coordinate y at various value of moving stream velocity 
U q are plotted in Figure 5 for two dimensional jet and in Figure 6 for 
axisymmetrical jet. For all the cases, the initial velocity profiles 
are rectangular. These velocity distributions approach to the form of 
linearized solution when U q is much larger than unity. 

When U Q « 0 , we compare the numerical solution with similar 

solutions of Eq, (14) or (15). To compare the axial velocity profile, 

we plot u/u m at constant { but different x In figure 7. For 

similar solution, u/u is a constant at constant £ but for numerical 

' m 

solution, u / u m depends on both £ and x . However as x increases, 
we find in Figure 7 that the numerical solution approaches to the similar 
solution. In Figures 8 and 9, the numerical solutions with two different 
initial velocity profiles are compared with the corresponding similar 
solutions for two dimensional and axisymmetrical jet respectively. In 
general, the similar solution agrees better with parabolic initial velocity 
distribution than with rectangular one. In all cases, the numerical solu- 
tion gives broader profile near the axis of the jet than the similar 
solution. This fact gives better agreement of numerical solution with 
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experimental data than the similar solution. 

(c) Transverse velocity profile v . When U Q ^ 0 , typical 

v - component velocity profiles at * - 0.1 and U Q - 1/2, 1, and 10 

for two dimensional jet and axisymmetrieal jet are shown In Figure IQ. 
When U Q >> 1, the v - velocity component Is indeed very small and can 
be neglected. 

When U Q * 0, Figures 11 and 12 show respectively the v - compo- 
nent at various x - station and for different initial velocity profiles 
for two dimensional and axisymmetrieal jet. The corresponding similar 
solution is also plotted In the same figure. For rectangular Initial 
profiles, the agreement between numerical and similar solutions becomes 
better as x q Increases. For parabolic Initial velocity distribution, 
the agreement is almost perfect for x > 0.2 . 

(d) Half width of jet. Figure 13 shows the spread of half 

width b * b /D where u/u « 1/2 as a function of axial distance x 

or x q for two dimensional and axisymmetrieal jet. Experimental data 

of reference 7 with proper correction as discussed above and also data of 

references 5 and 6 are also given in Figure 13. The experimental data 

show a broader jet than the theoretical prediction by about 10% for two 

dimensional jet and 2.5% for axisymmetrieal jet. The experimental line 

for rectangular initial profile of axisymmetrieal jet should not be com- 

pared directly with the numerical solution because there is a difference 

in the maximum jet velocity in these two cases because of the momentum 

loss in the experimental data. Should x * 0.3 be used as recommended 

s 
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In section (a), then the experimental data Is again broader than 
theoretical prediction as described In Reference 7. The spread of 
half width of linearized solution is also shown In figure 13. 


PART R. THREE DIMENSIONAL JET 
VII. FUNDAMENTAL EQUATIONS 


The non-dimensional fundamental equations for a steady three 

3 

dimensional laminar jet mixing of an incompressible fluid are as 
* 

follows ? 


2 2 

/ . \3u , 9u . 9u 9 u , 9 u 

( U o + U )_ + V + w ^ + - I 


(U + u )i!i + v iw + 3w 8fw + j^w 

lU o + U '3x + v 3y + w 3z 0y 2 + 0 2 


3u + + Jw n 

3x + 3y + 3 Z 0 


( 21 ) 


The non-dimensional quantities are the same as those of part A with 

* * 
w « / 

w » R “x and z » z /D . 

0 ft ' 

j 

The lntlal and boundary conditions for Eq. (21) are: 


*From a discussion of Professor Isaac Greber, the z-momentum equation 
should be replaced by a vortlclty equation. Since this change will not 
affect the results of the present report, we shall look into this problem 
in our future research program. 
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x * 0 i 
x > 0 ; 


where L la the aspect ratio of rectangular exit.. In this report, we 

shall consider the case U » 0 only. 

o 

VXXI. NUMERICAL METHOD 

The alternating direction implicit finite difference method described 
In reference 8 seems most suitable for the present problem. The basic 
principle Is the same as that described In section V. However, the 
progress In x - direction Is completed in two steps. First we treat 
one direction (e.g., z - direction) as unchanged and then we treat the 
other direction (e.g., y - direction) as unchanged. Therefore, In each 
step, the problem is essentially the same as two dimensional one. This 
method was proved to be unconditionally stable for linear three dimen- 
sional equation. 

Applying forward difference in x - direction and central diffe- 
rence in y and z directions as described in section y, the x-momentum 
equation in Eq, (21) can be written in finite different form as follows; 


u^ » 1 for y & 1 and z < 1 

u . * 0 for y > 1 , z > h 


v « 


w 


u 


3w 

3y 

3v 

9z 

3w 

dz 


0 

0 

3v 

dy 


along 
along 
* 0 at 


z ■ 0 


y * 0 


( 22 ) 


y * « , z «> . 
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C 1 R+ + C T“ + 2)u £-ii,m,n + ^ " 1 > u 




^ - r^l-lMhn - (2 - ¥> u £-l, m ,„ + Cl + C23a) 


,kv 


(f- + X) Up , + + 2)u» + (§£ - l) Uj , 

I £,m-l,n v h £,m,n ' 2 ' £,m+l,n 


kv 


k 

P 


A ’l 

2 1(1 _ f’ )u l-h, m,n+l ” (2 - T- )u £-% >m ,n + a + C23b) 


where h is the half spacing in the x - direction, and p, k are 
respectively the y- and z-spacing* The subscripts l, m, and n 
represent the position of grid point in x, y and z direction respectively, 
The finite difference forms for the z~momentum equation in Eq, (21) 
are the same as Eqs. (23a) and (23b) . They will not be written down. Be- 
cause we find that the numerical solutions for these complete non-linear 
three dimensional equations (23), etc. are highly unstable. We have to use 
some approximation described in the next section to obtain some useful 
results of the present problem in which the z - momentum equation may be 
neglected, 

XX. APPROXIMATE NUMERICAL SOLUTION 
Attempts were first made to solve Eq. (23) numerically with all 
u, y and w . Since the Initial values of v and w are unknown and they 
are assumed to be zero at the nozzle exit, we must calculate the values of 


y.iiiiniijiSi. 
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v and w at the first station by satisfying the z - momentum 
equation and equation of continuity, No stable numerical solutions 
of v and w have been obtained In this way and numerical integra- 
tion can not proceed downstream. We are still trying to improve our 
numerical method in order to find stable numerical solution including 

the cross flow term. However from our successful numerical solutions 

/ 

of two dimensional and axlsymmetrlcal jet reported In Part A, we find 
out that the Influence of cross flow terms is not too strong on the 
distlrbutlon of axial velocity distribution. If our main interest is 
the axial velocity distribution the following equation 



SK 


2 — 2 — 

a u , a u 

2 *** 2 

8y 


(24) 


together with the proper initial and boundary condition (22) would give 
reasonable good results. We use bar to denote without cross flow. Hence 
we solve equation (24) numerically for axial velocity distribution with 
the initial and boundary condition (22) . In order to check the accuracy 
of the approximation (24) , we calculate also equation (24) for the two 
limiting cases; (i) two dimensional case L - °° and (il) axisymmetrical 
case. Since we have the exact numerical solution for these two limiting 
cases in Part A, we may easily find out the accuracy of this approxima- 
tion (24) by comparing the numerical results of Eq„ (24) of these two 
limiting cases and the exact solutions with cross flow. We find the 
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accuracy is reasonable good. Since we expect that the general three 
dimensional results should lie between these two limiting cases, we 
may find a method of correction to improve the numerical results of 
three dimensional jet Of Eq. (24) which will be discussed later, 

X. RESULTS AND DISCUSSION 

(a) Numerical computation. In Part A, Ay « 0.05 and 
y/D * 20 were used in the computation and satisfactory ^results were 
obtained. In the computation of three dimensional .jet, due to the 
storage capacity of UNIVAC 1108, total grid point over a station can 
only have 70 x 70 . When Ax ~ Ay = 0.1 are used, we have z/D * y/D - 7. 
This means that the boundary is set at 7D away in the y- and z-direction. 
A check of using such a boundary and spacing in two dimensional and axi- 
symmetrical jets indicates poor result for x > 1 . The central (maximum) 
velocity falls too rapidly for x > 1 . However, for x < 1 , the agree- 
ment with previous solution is good. To improve this result, we increase 
Ay and Az to 0.2 and then the boundaries are at y/D * z/D = 14 . 

Thus, we get better boundary conditions but sacrifice the accuracy in 
finite difference method. A check for two dimensional and ax ^symmetric 
results shows that the center velocity is slightly lower by 0.01 to 
0.02 in magnitude and the velocity profile is lower by the same amount 
near the central portion but agrees well in the tail. Most important is 
that the whole flow field is in better agreement with the exact solution. 
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Therefore, In our three dimensional jet calculation, both values 
Ay * Az « 0.1 (for x $ 1 only) and Ay * Az * 0.2 were used. When 
we compare the numerical results for these two spacing values, it was 
found that the difference in center velocity is insignificant and the 
velocity profile shows a wider tail in z - direction when Ay * Az ■ 0.2, 
were used with x £ 0.3 and L > 1 but insignificant difference in the 
y - direction. Thus we feel that the spacing Ay » Az » 0.2 used in the 
numerical Integration is sufficient. The same program was applied to the 
three dimensional linearized equation and the numerical results agree 
satisfactorily with the analytic solution of reference 3. 

(b) Central velocity (maximum velocity at each x - station) . 

In Figure 14, numerical solution of the central velocity (without cross 

velocity) as a function of x - distance from nozzle exit are plotted for 

two dimensional, axlsymmetrical and three dimensional jet with aspect ratio 

1*1,2 and 4 without free stream. It is interesting to note that 

the character of these curves are similar to those obtained from linearized 
1 

equation . For L * 1 (also axlsymmetrical jet) and L = 2 , no significant 
difference between the numerical solution and linearized solution was found. 
For large L , numerical solution gives higher center velocity. In Figure 14, 
exact numerical solutions with cross velocity for two dimensional and axi- 
symmetrlcal jet are also shown. 

We may figure out the correction factor due to neglecting the cross 
flow for two dimensional and axisymmetric ' jets by comparing the corresponding 
approximate and exact numerical solutions. Furthermore, since the two dimen- 
sional and axlsymmetrical jets are two limiting cases of the three dimensional 
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jet, we may consider the correction factors of two dimensional and 
axisyimetrical jet as lower and upper bounds respectively for the 
correction factor of the three dimensional jet. The following correc- 
tion factor for three dimensional jet is suggested. 

Let u and u be the center velocity with cross flow 

mo m 00 

for axisymmetrical (subject o) and two dimensional (subscript °°) jet 

respectively. Then u , and u T as well as c , c , c r be 

mo mp° mL uo u°° uL 

the central velocity without cross flow and correction factors for axi- 
symmetrical, two dimensional and three dimensional (subscript L) jet 
with aspect ratio L respectively. We define 


and 
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mo 
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m 00 
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(26) 


Au = u — u 
mo m^ mo 





(27) 


Then 


Au 


mL 


uL Au 


mo 


Ac 4- c 
uo u 00 


(28) 




(1 + c 



(29) 


Figure 15 shows the correction factor c^ as a function of distance 
from nozzle exit. The corrected central velocity for L - 1, 2 and 4 are 

shown in Figure 14. 
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(c) Axial velocity profiles. The axial velocity profiles for 
three dimensional jet without cross flow are shown in Figure 16 for 
L ■ 1, 2 and 4 . In Figure 17, the errors introduced by neglecting 

the cross flow for two dimensional and axlsymmetrical jets are shown by 
plotting (u/u) / (u m /u ) where u and u denote axial velocity with and 
without cross flow respectively and subscript m denotes the value on the 
axis of the jet, i.e., u m is the central velocity. It should be noted 
that even though, the error at the tail is large but the effect is small 
because the magnitude of velocity at the tail is always very small. 

Similarly, we may find the correction factor for the general ve- 
locity profile from the two dimensional and axisymmetrical jet results in 
the same manner as that for center velocity. The difference in error due 
to cross flow for two dimensional and axisymmetrical jets 


A 


u 


u 


u 


u 


[fcr> / fc~)] - [fcr) / (~) 1 


u u u u 

o mo » m 00 


(30) 


is always small in the central portion but it may become large toward the 
tail where the magnitude of velocity is always small. Hence a correction 
curve based either on two dimensional or axisymmetrical jet should be suffi- 
ciently accurate for the three dimensional jet. 

The shape of three dimensional jet varies with x and their contours 
of constant velocity over a station are different for different velocity. For 
the purpose of correction, the local shape of the three dimensional jet plays 
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an important role. For simplicity, the contour of half-width (u^u^/2) 
is chosen as the representative shape of jet for all velocities over a 
given x-station. The contours of half-width for three dimensional jet 
without cross flow are shown in Figure 18, The following rules are 
suggested for the correction of three dimensional jet due to cross flow 
effects : 

(i) For plane perpendicular to y-z plane with largest half- 
width b^ m , behaves like an axi symmetrical jet. 

(ii) For other planes perpendicular to y-z plane half-width 

approaches to two dimensional jet linearly proportional to b^/b^ , 

we have 


(U L /U L ) 



/u mL ) 



(u 

' m 00 


/u ) 
m® 



(31) 


where A is given by Eq. (30). Since ( u mI / u mL ) is given in Eq. (29) and 
as well as all terms on the right-hand side of Eq. (31) are known, we 
may calculate the value u^ . We feel that this correction is sufficient 
accurate even though the correction scheme is somewhat arbitrary. 

(d) Half-width of the three dimensional jet. The contours of 
half-width of three dimensional jet with L - 1, 2 and 4 without cross 
flow are shown in Figure 18. The correction factors due to cross flow are 
shown In Figure 15 which were obtained in a similar manner as other correc- 
tion factors from the exact solutions of two dimensional and axisymmetrical 
jet. The difference of these correction factors for the two limiting cases 


is small. 
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XX. CONCLUSIONS AND GENERAL REMARKS 

From our numerical solutions, the following conclusions may 
be drawn: 

(1) Numerical solutions for two dimensional and axlsymmetrical 
free jets give better description of the flow field than similarity so- 
lutions as compared with the existing experimental data. The boundary 
layer equation is a good approximation for the free jet problem provided 
that the Reynolds number at nozzle exit is not too low* 

(2) For two dimensional and axlsymmetrical jets, parabolic 
initial velocity distribution gives better agreement between numerical 
and similarity solutions. The approach of numerical solution at large 

x to similarity solution is demonstrated clearly in both central velocity 
and velocity profile. 

(3) When the jet velocity deviates slightly from the uniform 
free stream velocity, the numerical solutions agree with the results of 
linearized theory. For instance, when the jet velocity is 1/10 greater 
than the free stream value, the difference between numerical and linearized 
solutions is of the order of 0.001 to 0.002 of the free stream velocity. 
With free stream, the magnitude of the cross flow is reduced significantly. 

(4) The axial velocity distributions of the three dimensional jet 
issuing from a rectangular nozzle of various aspect ratio without free 
stream have been calculated numerically by neglecting the cross velocity 
(v * w 3 0) • The correction factors due to cross flow have been obtained 
from the exact solutions with cross flow of two dimensional and axlsymmetrical 


jet. 
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(5) The central velocity given by non-linear equations is 
higher than the corresponding linearized equation and the spread of 
non-linear jet is much greater than the linearized case. 

(6) The effects of neglecting the cross flow are to lower the 

i 

center velocity and to increase the velocity in the tail. 

(7) Although no experimental data available for the three 
dimensional jet, from the results of the two dimensional and the 
axlsymmetrical jet, we feel that our three dimensional jet numerical 
solutions are sufficiently accurate. 

(8) Further development of numerical scheme to include the cross 
flow for three dimensional jet is still needed, because in our approximate 
numerical solution, no information about cross velocity components for 
three dimensional jet is given. 
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APPENDIX I 


In the Implicit finite difference method described In Part A for 
two dimensional and round jet and Part B for three dimensional jet, u- 
component In the new station must be obtained by solving a set of N-l 
simultaneous algebraic equations as follow; 


A rVl + Vr + c r u r+l " D r 
A N-X U N-2 + " D N-1 


r «« 1 

1 < r < N-2 
r - N-l 


where N Is the total number of grid points In one direction, A, B, C 
and D are coefficients. Treat A, B, C and D as vectors of N-l component 
and let a, 3 and y be the new vectors of N-l component and 





r-1 


Q "I 

3 e a c 
r r r 


■ °i D i 


-i 


Y * a (D - y -A ) 
? r r N r 1 r-1 r 


2 < r * N-l 
1 £ r $ N-l 


2 ^ r $ N-l 
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it can be shown that 


Vi " y n~i 



3 u .. 
r r+1 


1 $ r 


Thus a , 8, y are calculated in order of increasing 


$ N-2 


and u 


is obtained in order of decreasing r . 


For proof, see Reference 4. 
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A computer program for two dimensional jet with rectangular 
initial velocity distribution is given below; 


C Two Dimensional Jet, Numerical Calculation (Implicit Method) 

Dimension U(800), UP(800), B(800), V(800) , AM(800) , AN (800) , 
AU(800), ID(8G0), AHI(10), LI(10), A(800) , C(800) 

C STEP INITIAL VEL. DISTRIBUTION AT THE EXIT 

9900 RE AD (5, 1100) NY, NYE, LL* KY, KYI, KY2, KY3, LY, P, UO, XST 

WRITE (6,1100) NY, NYE, LL, KY, KYI, KY2, KY3, LY, P, UO, XST 

1100 FORMAT (815, 3F10,5) 

READ (5,3100) (LI(N), N-l, LL) 

WRITE (6, 3100) (LI(N) , N-l, LL) 

3100 FORMAT (1015) 

READ(5,3200) (AHI(N) , N-l, LL) 

WRITE (6, 3200) (AHI(N), N-l, LL) 

3200 FORMAT (3F10, 6, 7F6.5) 

NY1 - NY + 1 

NYR - NX - 1 

KD * KY - 1 

DO 51 I * 1, NY1 

IF(I-NYE) 501, 501, 503 

501 UP (I) - 1 

GO TO 51 

503 UP(1) » 0. 

51 CONTINUE 

WRITE (6, 1009) 

1009 FORMAT (30H INITIAL VELOCITY DISTRIBUTION,/) 

WRITE (6, 1006) 

WRITE (6 ,1000) (UP (I) , I - 1, NY1, LY) 

WRITE (6, 1001) 

X - 0 

UP (NYE) * .5 
P2 - P**2 
KYP » 10 

DO 9800 KL » 1, LL 
L = LI(KL) 

AH * AHI(KL) 

DO 9000 KKK - 1,L 
X « X + AH 
IF(UO) 81, 82, 81 
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82 

81 

83 

201 

202 

203 

17 


15 

306 

61 

305 

14 

9100 

63 

29 


XPM * 0 
GO TO 83 

XPM «* X/UO 

00 9200 K * 1, KYP 
DO 9100 KA * 1, KYI 
DO 17 H * 1, NY1 
IP(KA - 1) 201, 201, 202 

UCO-UP(M) + UO 
GO TO 203 

UC0®U(M) 4* UO 

B(M)H?2*UCO/AH+2, 

D (M) v UCO *P 2 *UP (M) / AH #1 
AM(1)*B(1) 

AN(1)*~2. /AM(1) 

AU (1) «D (.1) / AM( 1) 

DO 15 M-2.NY nr> 

AA*-(.5*P*V(M) + 1) Jr 
CA«(, 5*P*V(M) - 1) 
AM(M)* b B(M) - AA*AN (M~l) 
AN(M)*CA/AM(M) 


AU(M)»(D(M) - AU(M-1)*AA)/AM(M) 
XF(KA - KYI) 305, 306, 306 

DO 61 M ~ 1, NY1 

B(M) * U(M) 

U(NY) » AO (NY) 

DO 14 M * 1, NYR 
MM - NY - M 

U (MM) *AU (MM) -AN (MM) *U (MMf 1) 

CONTINUE 
DO 63 M » 1, NY 

U(M) * , 5*(B(M) + U(M>) 

V(D - o. #4 

SUM * (U(l) - UP (1)) /Ail f 
DO 29 M * 2, NY1 ,, 

DDT«(U(M) - UP (M) ) /Air 
SUM'SUM + 2.*DDT n, 

V(M) * -(SUM - DDT)*.5*P 
KYI « KY2 

IF (KYP - 1) 9200, 9200, 309 
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309 

308 

67 

9200 

68 
1016 

601 

65 

402 

9000 

401 

9800 

9500 

1000 

1001 

1006 

1007 


IF(K-KD) 9200, 308, 9200 

DO 67 M-l, NY1 
C(M) * V(M) 

A(M) - U(M) 

CONTINUE 

DO 68 M * 1, NY1 
V(M) « ,5*(C(M) + V(M) ) 

U(M) - .5*(U(M) + A(M)) 

WRITE (6, 1016) X, XPM 

FORMAT ( / , 30H. AVERAGE VALUE FOR U, V AT X*,F9 .5 ,5X,4H X 
WRITE (6, 1006) 

WRITE (6 ,1000) (U(I) ,I«1,NY1,LY) 

WRITE (6, 1007) 

WRITE (6 , 1000) (V (I) ,1*1 ,NY1 ,LY) 


DO 65 M-l, NY1 
UP (M) * U(M) 

IF(X-XST) 9000, 402, 402 

KYP * KY 

CONTINUE 

IF(KL-l) 9800, 401, 401 

KY1-KY3 

KY2-KY3 


CONTINUE 
WRITE (6, 9500) 

FORMAT (1H1) 

GO TO 9900 

FORMAT (10F12 * 5) 

FORMAT (/ /) 

FORMAT (12R U COMPONENT, /) 
FORMAT ( / ,12H V COMPONENT,/) 


-,F9.5,/) 


END 
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where NY » total number of grid point in y direction minus 1. 

NYE » grid point at edge of exit 

LL * number of sets of equal spacing stations to be calculated 

KY « number of iteration for v at downstream X > XST 

KYI, KY2, KY3 ■ number of iteration for u at first section* first 

set of stations and following stations respectively 

LY ® grid point in interval for printing 

p ~ spacing in y direction 

U4> 53 free stream velocity (U q /Uj) 

„ * 

XST ■ distance from exit, for X < XST iteration in v is 10; 

* 

for X < XST iteration in v is given by KY. 

LI(LL) = number of stations to be calculated in each set 
AHI(LL) » spacing between stations in each set 

Transform from two dimensional case to axisymmetrical case can easily be 
done by changing the statements with # sign as follow: 

#1 By adding a statement immediately next to it 

B(l) « B{1) + 2 

n AA « — ( . 5*P*(V(M) - 1/((M-1)*P)) + 1) 

#3 CA = ._^V'V(M) - 1/ ( (M-l) *P) ) - 1 

#4 MULTIPLYING BY 0.5 ON THE R.H.S. 

#5 MULTIPLYING BY (M-l) ON THE R.H.S. 

#6 DIVIDING BY (M-l) ON THE R.H.S. 









001 .01 .10 1.0 10.0 

X = vX*/U,*D 2 

FIG. 3. CENTRAL VELOCITY AS A FUNCTION OF DISTANCE FOR VARIOUS INITIAL 
VELOCITY DISTRIBUTIONS AT EXIT, TWO DIMENSIONAL JET WITH AM- 
BIENT FLUID STILL. 
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ROUND JET WITH MOVING STREAM 



FIG. 6. VELOCITY PROFILES FOR AXISYMMETRICAL JET WITH MOVING STREAM. 




(a) TWO DIMENSIONAL JET (b) AXISYMMETRICAL JET 

FIG. 7. COMPARISON OF u/u m BETWEEN NUMERICAL AND SIMILAR SOLUTION AS A 
FUNCTION OF DISTANCE FROM EXIT 
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NUMERICAL SOLUTION 



(A). RECTANGULAR INITIAL VELOCITY DISTRIBUTION. 

FIG. 9. COMPARISON OF VELOCITY PROFILES BETWEEN NUMERICAL AND 
SIMILAR SOLUTION AT CONSTANT MOMENTUM FOR AXISYMMETRICAL 
JET WITH AMBIENT FLUID STILL. 
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FIG. 10. VERTICAL VELOCITY DISTRIBUTION AT X 0 =O.I 
FOR JETS WITH MOVING STREAM, RECTANGU- 
LAR INITIAL VELOCITY DISTRIBUTION. 
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FIG. II. VERTICAL VELOCITY DISTRIBUTION FOR TWO 

DIMENSIONAL JET WITH AMBIENT FLUID STILL. 
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FIG. 12. VERTICAL VELOCITY DISTRIBUTION FOR 

AXISYMMETRICAL JET WITH AMBIENT FLUID 
STILL. 


±s 

a: o 
ujo: 
cl u. 
x Vt 

£q 

t uj 
>2 

25 
O CE 

3 § 

o Q 
to >- 
□ 

£T -J 

< <ro 

= Zu: 
5 UJ uj 
to S £E 




cn 

I- 

UJ 

o 

x 

H 

Q 


Ll. 

< 

X 

Ll. 

O 

Q 

< 

UJ 

X 

X 

if) 

UJ 

X 

I— 

K) 

CD 

iL. 




AXISY. JET 



FIG. 15. THE CORRECTION FACTOR Cu L AND HALFWIDTH RATIO b/b AS 
FUNCTION OF DOWNSTREAM DISTANCE FROM EXIT. 




FIG. 16. VELOCITY PROFILE FOR THREE DIMENSIONAL JET. 




n/ n) / (n/ 


5 


6 


THREE DIMEN- 
FLOW. 




